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ABSTRACT 

q_j Fundamental parameters and the carbon, nitrogen and oxygen abundances are 

determined for 22 B-type stars with distances d ^ 600 pc and slow rotation (v sin i ^ 
66 km s _1 ). The stars arc selected according to their effective temperatures T e g and 
surface gravities logg, namely: T e g is between 15300 and 24100 K and log<? is mostly 
greater than 3.75; therefore, stars with medium masses of 5-11 M Q are selected. The- 
ory predicts for the stars with such parameters that the C, N and O abundances in 
their atmospheres should correspond to their initial values. Non-LTE analysis of C II, 
N II and O II lines is implemented. The following mean C, N and O abundances are 
obtained: log e(C) = 8.31±0.13, log e(N) = 7.80±0.12 and log e(O) = 8.73±0.13. These 
values are in very good agreement with recent data on the C, N and O abundances for 
nearby B stars from other authors; it is important that different techniques are applied 
by us and other authors. When excluding for the stars HR 1810 and HR 2938, which 
can be mixed, we obtain the following mean abundances for the remaining 20 stars: 
loge(C) = 8.33±0.11, loge(N) = 7.78±0.09 and loge(O) = 8.72±0.12; these values 
arc in excellent agreement with a present-day Cosmic Abundance Standard (CAS) of 
Nieva & Przybilla. 

The derived mean N and O abundances in unevolved B stars are very close to 
the solar photospheric abundances, as well as to the protosolar ones. However, the 
mean C abundance is somewhat lower than the solar one; this small but stable carbon 
deficiency is confirmed by other authors. One may suggest two possibilities to explain 
the observed C deficiency. First, current non-LTE computations of C II lines arc still 
partially inadequate. In this case the C deficiency is invalid, so one may conclude that 
the Sun and the local unevolved B stars have the same metallicity. This would mean 
that during the Sun's life (i.e., for the past 4.5 • 10 9 yr) the metallicity of the solar 
neighbourhood has not markedly changed; so, an intensive enrichment of the solar 
neighbourhood by metals occurred before the Sun's birth. Second, the C deficiency in 
the local B stars is valid; it is supposed that the Sun can migrate during its life from 
inner parts of the Galactic disk where it has born, so its observed chemical composition 
can differ from the composition of young stars in its present neighbourhood. 

Key words: stars: abundances - stars: evolution - stars: early-type 



1 INTRODUCTION 

Some recent works show that the metallicity of young stars 
in the solar neighborhood is very close to the Sun's metallic- 
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ity. In particular, it was shown by Lyubimkov et al. (2010) 
that the mean iron abundance in the atmosphere of 48 A- 
, F- and G-type supergiants with distances d < 700 pc is 
log e(Fe) = 7.48±0.09, a value that coincides with the solar 
abundance loge (Fe) = 7.50±0.04 (Asplund et al. 2009). 
Earlier it was found for local early B-type stars, the progen- 
itors of AFG supergiants, that their mean magnesium abun- 
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dance is loge(Mg) = 7.59±0.15 (Lyubimkov et al. 2005), 
that is equal to the solar abundance log £0(Mg) = 7.60±0.04 
(Asplund et al. 2009). One may mention also the Fe abun- 
dance obtained by Luck, Kovtyukh & Andrievsky (2006) 
for Cepheids in the solar vicinity, as well as the Fe and Mg 
abundances found by Fuhrmann (2004) for nearby F, G and 
K dwarfs and subgiants of the thin Galactic disk. Moreover, 
it was found for B-type stars in the Orion OBI association 
that their S and O abundances are very close to the solar 
ones; see Daflon et al. (2009) and Simon-Di'az (2010), respec- 
tively. Nieva & Simon-Di'az (2011) showed that the B stars 
in the Orion OBI association show good agreement with the 
Sun in their C, N, O. Mg, Si, and Fe abundances. 

Just three years ago, the problem of metallicity of young 
stars seemed far from simple. It was shown by several au- 
thors that the metal abundances in nearby B-type stars were 
somewhat less than in the Sun (see, e.g., Morel's 2009 re- 
view). In particular, the derived C, N and O abundances of 
B-type stars tended to be less than in the Sun. Two possi- 
ble explanations were considered, namely: (i) there is a real 
slight deficiency of C, N and O in the stellar atmospheres, 
e.g., connected with stellar evolution or (ii) the C, N and 
O abundances were systematically underestimated in some 
studies. Then in 2003-2008, studies were reported where the 
N and O abundances in B stars were close to the revised N 
and O solar abundances based on a three-dimensional model 
of the Sun's atmosphere; the C abundance continued to show 
a small deficiency (see, e.g., Herrero 2003; Herrero & Lennon 
2004 and Przybilla 2008). In order to investigate the possi- 
ble metallicity difference between the Sun and young stars, 
we undertook an accurate analysis of the C, N and O abun- 
dances for a rather large sample of nearby early and medium 
B-type stars in the main sequence (MS) evolutionary phase. 
(Note that up-to-date results for C, N and O appeared dur- 
ing our analyses will be discussed in Section 7). 

The C, N and O abundances in early and medium B- 
type stars are determined here from C II, N II and O II spec- 
tral lines. These lines are rather sensitive to two fundamental 
parameters of the stars, namely their effective temperatures 
T c ff and surface gravities log g. Recently, we showed (Lyu- 
bimkov et al. 2009, 2010) that a significant improvement in 
the accuracy of the log g values can be obtained through ap- 
plication of van Leeuwen's (2007) new reduction of the Hip- 
parcos parallaxes. Using these new parallaxes in an analysis 
of A, F and G supergiants, we obtained logg with very high 
accuracy for stars with distances d < 700 pc. Now, we apply 
this method in the present study of B-type MS stars. 

As far as the effective temperature T c ff is concerned, 
there is some difference in various T a g scales for B stars 
that is especially significant for the hottest B stars with 
T c ff ^ 25000 K. Therefore, we consider relatively cool B 
stars with T c ff ^ 24000 K, for which the Teff estimates are 
more reliable. Moreover, we consider mostly those B stars 
with \ogg ^ 3.75, i.e., stars rather far from the termination 
of the MS phase, in order to exclude possible evolutionary 
changes in the C, N and O abundances. In other words, we 
try to determine the initial abundances of these elements in 
such young stars. 

Carbon, nitrogen and oxygen participate in the CNO- 
cycle, which is the main source of energy of the B-type MS 
stars. The C and N abundances change significantly in stel- 
lar interiors during the MS phase (the O abundance shows 



smaller alterations). It has been predicted that the C, N and 
O abundances can change markedly in surface layers too, if 
rotationally-induced mixing in the MS phase occurs (Heger 
& Langer, 2000; Frischknecht et al. 2010). The changes in 
the surface C, N and O abundances by the MS phase ter- 
mination depend strongly on the mass M and the initial 
rotational velocity vo of a star: the greater M and vo the 
stronger the predicted final alterations. Since our goal is to 
obtain the initial abundances, we consider the B-type stars 
with medium masses (5-11 M@) that are not close to the 
MS phase termination. Their observed rotational velocities 
are rather small, too (vsini ^ 66 km s _1 ) which implies 
that the majority of the selected stars are slow rotators but 
an admixture of rapid rotators observed at high angles of 
inclination may be present. 



2 SELECTION OF STARS 

We have studied previously a large sample of early and 
medium B-type stars in the MS phase; results have been 
published in four papers (Lyubimkov et al. 2000, 2002, 2004 
and 2005; hereinafter Papers I, II, III and IV, respectively). 
In particular, in Paper III we determined from He I lines 
the helium abundance He/H and the projected rotational 
velocity v sin i for 102 B stars. The v sin i values were found 
to range from 5 to 280 km s _1 ; such a large variation allowed 
us to analyze a relation between the helium abundances and 
rotational velocities. The used He I lines are strong in spec- 
tra of early and medium B stars, so their equivalent widths 
can be measured without problems even for the stars with 
rather high vsini. On the contrary, the C II, N II and O II 
lines used in the C, N and O abundance analysis are signif- 
icantly weaker, so measurements of their equivalent widths 
for vsini > 100 km s _1 are difficult, especially for stars with 
effective temperatures T a a < 17000 K. As a result, we found 
that only 54 stars with v sini < 100 km s _1 from our original 
list are suitable for the C, N and O abundance analysis. 

There are other limitations as well. One concerns the 
effective temperatures T c ff of B stars. It is known that the 
T e ff scales of various authors for B stars are different, and 
the difference increases with T c g , so the differences are espe- 
cially large for the hottest B stars with T e ff ^ 25000 K (see, 
e.g., a comparison of some T e ff scales in Daflon, Cunha & 
Butler 2004). We determined in Paper II the T e ff and logp 
values for more than 100 B stars and confirmed that uncer- 
tainties in the T e ff are especially great for T c ff ^ 25000 K. 
Therefore, we selected from our original list only stars with 
Teff ^ 24000 K. The next limitation regards another funda- 
mental parameter, the surface gravity log g. In order to ex- 
clude as far as possible evolutionary effects on the observed 
C, N and O abundances, we consider stars which are not 
close to the MS phase termination. Therefore, the selected 
stars have mostly logp > 3.75. 

One more limitation regards the stellar parallaxes 7r, 
which are used by us as the principal indicator of the 
surface gravity logg. Using van Leeuwen's (2007) catalog, 
we selected stars with ir ^ 1.67 mas, i.e., with distances 
d ^ 600 pc (mas = milliarcsecond) . Selection of stars with 
smaller n values would reduce markedly the accuracy of the 
derived surface gravities log p. 

The final list of B stars selected is presented in Table 1; 
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Table 1. Basic parameters of programme B stars 



HR 


HD 




TV, 

mas 




d, 
pc 


V s'mi. 
km s — 


l 


T c ff 
K 




logg 

cgs 




M/Mq 


38 


829 


1 


81±0 


.42 


552 ±98 


11.5±1 


.5 


19220±510 


3 


83±0. 


15 


7.6 


1617 


32249 


4 


,42±0 


.23 


226 ±11 


45.0±1 





18570±20 


3 


,87±0. 


04 


7.1 


1640 


32612 


1 


76±0 


.41 


568±102 


53.5±1 


.5 


19750±230 


3 


,81±0. 


15 


8.0 


1781 


35299 


3 


72±0 


.32 


269 ±21 


5.0 ±2. 


5 


23700±410 


1 


,27±0. 


07 


8.9 


1810 


35708 


5 


,20±0 


.21 


192 ± 7 


24.0±2 


.0 


21040±520 


4 


11±0. 


07 


7.7 


1820 


35912 


2 


,53±0 


.55 


395 ±67 


13.5±2 


.0 


19370±380 


4. 


,05±0. 


19 


7.0 


1923 


37356 


2 


,03±0 


.47 


493 ±88 


17.5±3 


.0 


21400±630 


3 


,78±0. 


16 


9.5 


1933 


37481 


2 


44±0. 


39 


410 ±55 


66.0±2 


.5 


24090±850 


4. 


,14±0. 


13 


9.7 


2058 


39777 


2 


80±0. 


53 


357 ±55 


21.5±2 


.0 


21720±370 


4 


24±0. 


16 


7.7 


2205 


42690 


2 


73±0. 


.28 


366 ±34 


9.8 ±3. 


5 


19390±310 


3 


,74±0. 


01 


8.1 


2344 


45546 


3 


00±0 


.24 


333 ±25 


61.0±2 




19420±570 


3 


,76±0. 


07 


8.0 


2756 


56342 


5 


.19±0. 


22 


193 ± 8 


26.0±1 




16590±390 


4. 


06±0. 


04 


5,1 


2824 


58325 


1 


67±0 


36 


599±101 


12.5±2 


.5 


20210±450 


3 


86±0. 


16 


8.1 


2928 


61068 


1 


93±0 


.28 


518 ±64 


26.0±1 




23650±850 


3 


,83±0. 


11 


11.2 


3023 


63271 


1 


99±0 


33 


503 ±69 


42.0±1 




22400±200 


3 


86±0. 


12 


9.7 


7426 


184171 


3 


80±0 


16 


263 ±11 


29.5±2 


.0 


16540±850 


3 


60±0. 


07 


6.8 


7862 


196035 


3 


56±0. 


46 


281 ±32 


34.0±1 


.5 


17750±550 


4. 


19±0. 


09 


5.7 


7996 


198820 


1 


87±0 


33 


535 ±78 


35.0±1 


.5 


15980±400 


3 


60±0. 


12 


6.5 


8385 


209008 


3 


07±0 


42 


326 ±38 


19.0±3 


.0 


15360±240 


3 


,76±0. 


10 


5.6 


8549 


212883 


2 


17±0 


.63 


461 ±95 


8.0 ±3. 





20400±680 


3 


,97±0. 


21 


7.8 


8768 


217811 


3 


07±0. 


58 


326 ±50 


8.0 ±4. 





18090±730 


3 


94±0. 


16 


6.5 


9005 


223128 


2 


19±0 


.27 


457 ±49 


15.5±2 


.5 


22340±870 


3 


79±0. 


11 


10.5 



it contains 22 objects in all. We present there for each star 
its HR and HD numbers, parallax n with its error from 
van Leeuwen's (2007) catalog and corresponding distance 
d = 1/n, where tt is in arcseconds. The presented errors 
in distances d are evaluated with the formula given by van 
Leeuwen's (2007, p. 86). The projected rotational velocities 
vsini with their uncertainties from Paper III are provided, 
too; note that the v sin i values have been found in Paper III 
from profiles of six He I lines. Other parameters in Table 1, 
namely the effective temperature T e g, surface gravity \ogg 
and mass M, will be determined below. 

The masses M of the selected stars are between 5 
and 11 Mq, and their projected rotational velocities are 
vsini ^ 66 km s~ . Therefore, we may hope that possible 
evolutionary changes to the surface N and O abundances 
due to the rotationally-induced mixing in the MS phase are 
minimized. However, we need to be confident that the de- 
rived N and O abundances of these stars are really close 
to their initial abundances; in other words, they have not 
been altered by mixing. In this connection, the stars with the 
lowest surface gravities logg are of special interest, because 
they are closer to the MS phase termination than other pro- 
gramme stars and, therefore, may show somewhat altered 
C, N and O abundances. There are seven stars in Table 1 
with log g < 3.80, specifically with logg = 3.60-3.79. Their 
masses are M = 5.6-10.5 Mq; the projected rotational veloc- 
ities are vsini = 9-35 km s _1 for six stars and usini = 61 
km s _1 for HR 2344. It follows from Frischknecht et al.'s 
(2010) computations that for the model with M = 12 Mq 
and initial rotational velocity vq — 62 km s _1 the surface 
value log e(C) decreases only by 0.02 dex and log e(N) in- 
creases only by 0.05 dex by the MS end. Therefore, for the 
lower masses and rotational velocities, like that for 7 stars in 
question, the expected changes in log e(C) and log e(N) are 
still smaller, whereas the changes in log e(O) are impercepti- 
ble (according to the theory, oxygen shows smaller changes 
than carbon and nitrogen. So, the derived loge(C), loge(N) 



and log e(O) values for all programme stars are anticipated 
to be very close to the initial abundances. 



3 OBSERVATIONS AND EQUIVALENT 
WIDTHS 

We used here the same high-resolution spectra as in Papers 
I-IV. The spectra were acquired in 1996-1998 at two obser- 
vatories, namely the McDonald Observatory (McDO) of the 
University of Texas and the Crimean Astrophysical Obser- 
vatory (CrAO). At the McDO the 2.7-m telescope and the 
Tull coude echelle spectrometer was employed (Tull et al. 
1995). The resolving power was R = 60000 and the typical 
signal-to-noise ratio was between 100 and 300. At the CrAO 
we observed on the 2.6-m telescope with the coude spectro- 
graph. In this case we had R = 30000 and a signal-to-noise 
ratio between 50 and 200. A more detailed description of 
the observations and reductions of spectra was presented in 
Paper I. It should be noted that spectral intervals of about 
70 A centered on He I lines were recorded at the CrAO. 
Since these intervals contain few C II, N II and O II lines, 
we used mostly the McDO echelle spectra for selection of 
the lines and measurements of their equivalent widths. 

When measuring the equivalent widths W , we found 
that rarely could the C II, N II and O II lines be accurately 
fitted by Gaussian profiles. Therefore, in most cases the 
W values were measured by direct integration. We checked 
for possible blending of the lines using the database VALD 
(Kupka et al. 1999, Heiter et al. 2008) and predicted W val- 
ues for models based on parameters T e g and log g from Paper 
II. Only those lines were selected where the estimated blend- 
ing was markedly less than 10 per cent. It should be noted 
that the number of measured lines is different for differ- 
ent stars; this number depends on the effective temperature 
T B ff and the observed rotational velocity vsini. The lines in 
question are especially weak for stars with T e g < 17000 K, 
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Figure 1. Comparison of our equivalent widths (in mA) of N II 
and O II lines with data from other works, namely: Gies & Lam- 
bert (1992) - squares, Cunha & Lambert (1992, 1994) - filled 
circles, and Kilian & Nissen (1989) - crosses. Straight line corre- 
sponds to the equal W values. 

so in this case the W values are measurable only for stars 
with low v sin i. For some stars it was possible to obtain W 
only for a few lines including the relatively strong lines C II 
4267 A and N II 3995 A. We found that use of these lines 
in the C and N abundance determination is questionable, so 
they were excluded from the analysis (see below). 

It is interesting to compare the measured equivalent 
widths W with other measurements. In Fig. 1 we compare 
our W values with the data from three sources, namely Gies 
& Lambert (1992), Cunha & Lambert (1992, 1994) and Kil- 
ian & Nissen (1989), where there are 3, 4 and 4 stars in 
common, respectively. One may see that there is good agree- 
ment for lines with W < 120 mA, and only the few lines 
with W > 120 mA show a small systematic discrepancy. So, 
one may conclude that there is, in general, good agreement 
between our and previous W measurements. 



4 NON-LTE COMPUTATIONS 

4.1 Necessity for a non-LTE approach 

Departures from local thermodynamic equilibrium (LTE) 
have an influence on computations of C II, N II and O II 
lines for early B stars. Using both LTE and non-LTE ap- 
proach, Gies & Lambert (1992) found that differences in 
the N and O abundances between these two cases are rather 
small on average, but for some stars they exceed 0.2 dex. An 
appreciable difference due to the non-LTE effects can appear 
for the derived microturbulent parameter Vt- According to 
Korotin et al. (1999b), non-LTE corrections to the oxygen 
abundance in early B stars can attain 0.2-0.3 dex. 

In order to study the non-LTE effect on C II, N II and 
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Figure 2. Comparison of computed equivalent widths of the N II 
4041.3 A line for non-LTE (solid curve) and LTE (dashed curve). 
Computations are implemented for (a) logg = 4.0 and (b) log g = 
3.5. 



O II lines in more detail, we performed both LTE and non- 
LTE calculations for all the lines used in our analysis - see 
the next section for details. We show as an example in Fig. 2 
and Fig. 3 results of the calculations for two lines, namely 
N II 4041.3 A and O II 4414.9 A (the solar N and O abun- 
dances are adopted). Equivalent widths W of the lines are 
displayed as a function of the effective temperature T e ff for 
(a) logg = 4.0 and (b) logg = 3.5. For these lines, as well 
as for other used lines, the W values in the non-LTE case 
are greater than in the LTE one. The difference is small for 
stars with T B ff < 18000 K but increases for hotter stars. 
Fig. 2 shows that in the W maximum region of the line 
N II 4041.3 A, i.e. at T cff = 25000 K (a) or 24000 K (b), 
the difference is 18 per cent for logp = 4.0 and 30 per 
cent for loggr = 3.5. According to Fig. 3, the discrepancy 
between non-LTE and LTE for the O II 4414.9 A line is 
more significant. In particular, in the W maximum region 
at T cB = 28000 K (a) or 26000 K (b) the difference is about 
40 per cent for log g = 4.0 and 55 per cent for log g — 3.5. So, 
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Figure 3. Comparison of computed equivalent widths of the O II 
4414.9 A line for non-LTE (solid curve) and LTE (dashed curve). 
Computations are implemented for (a) logg = 4.0 and (b) logg = 
3.5. 

it follows from our computations that departures from LTE 
for B stars with T e g > 20000 K are substantial. Therefore, 
we use a non-LTE approach for all lines and all stars in the 
present work. 

Our analysis showed that the non-LTE influence on de- 
rived abundances leads to a significant decrease of the micro- 
turbulent parameter Vt, especially for relatively hot B stars 
with T cfl > 20000 K. The non-LTE V t values for such stars 
can decrease by a factor of about two as compared with the 
LTE ones. Therefore, we confirmed Gies & Lambert's (1992) 
conclusion about the strong dependence of the Vt derivation 
on the LTE or non-LTE approach. 

4.2 Non-LTE computations and selection of lines 

For the non-LTE computations of C II, N II and O II lines, 
we applied the code MULTI created by Carlsson (1986) and 
updated by S.A.Korotin (see Korotin et al. 1998, 1999 a,b). 
This code, as well as all input data including the N and 
O model atoms was kindly provided to us by S.A. Korotin 



(the carbon model atom is taken from Sigut, 1996). Note 
that we used this same code in a non-LTE analysis of the 
nitrogen abundance from N I lines for A- and F-type super- 
giants (Lyubimkov et al. 2011). Atomic data for C II, N II 
and O II lines are taken from the database VALD (Kupka et 
al. 1999; Heiter et al. 2008). The list of the lines used is pre- 
sented in Table 2 including their lower excitation potentials 
Ei and the oscillator strengths log gf. In the computations 
we used, as in Paper III and IV, the plane-parallel blanketed 
LTE model atmospheres, which were computed by us with 
the code ATLAS9 (Kurucz 1993) for the corresponding pa- 
rameters T c ft and logp of programme stars. Our approach 
is what has been termed a hybrid method, i.e., LTE model 
atmospheres are combined with non-LTE line formation. 

As mentioned above, 22 B stars with effective temper- 
atures T fr ^ 24000 K were selected for the present work. 
However, we plan to analyze later the C, N and O abun- 
dances for hotter B stars from our original list with T e g up 
to 30000 K. Therefore, from the outset, we performed non- 
LTE computations of C II, N II and O II lines for B stars 
with effective temperatures T e g up to 30000 K and rotational 
velocities vs'mi < 100 km s _1 . As a result, we considered 54 
objects in all instead of 22. Their parameters T B ff and logp 
were taken from Paper II as a first approximation. 

We compared for each of the 54 stars the non-LTE C, 
N and O abundances derived from individual C II, N II and 
O II lines with the mean C, N and O abundances obtained 
from other lines. We found that some lines display a sys- 
tematic discrepancy. In particular, the strongest C II line 
in the visible spectral region, namely the line C II 4267 A, 
showed a systematic difference in log e(C) for hotter stars of 
our sample: for stars with T e s between 20000 and 24000 K 
the log e(C) underestimation from 4267 A is about 0.2-0.7 
dex. A tendency of 4267 A to lower the loge(C) values for 
early B stars has been noted earlier. Meantime, a few at- 
tempts have been made to update the non-LTE calculations 
of this line by using a more detailed model atom. In par- 
ticular, Sigut 's (1996) computations improved markedly the 
agreement with observed equivalent widths of the line. We 
applied Sigut 's model atom in our non-LTE calculations of 
C II lines. Nevertheless, it was impossible to eliminate the 
significant discrepancy between 4267 A and other C II lines. 

Furthermore, a systematic difference was found for the 
strongest N II line in the visible spectral region: this line 
at N II 3995 A is frequently used in analyses of B and 
O stars. In Fig. 4 the difference between the log e(N) val- 
ues determined from the N II 3995 A line and the mean N 
abundances (without 3995 A) is shown as a function of T c g. 
The underestimate is obviously correlated with T e ff with the 
largest underestimation (up to 0.2-0.4 dex) occurring for the 
coolest B stars of the sample. If included in the collection 
of N II lines, this line has a significant influence on the de- 
termination of the microturbulent parameter Vt- It is clear 
that the 3995 A line should be excluded from the analysis 
of our programme stars with T e g ^ 24000 K. Nevertheless, 
as Fig. 4 shows, for the stars with T c s > 26000 K this line 
displays no systematic discrepancy; therefore, for such hot 
stars use of 3995 A in the N abundance analysis is likely to 
be permissible. 

Thus, we excluded from further analyses the C II 4267 A 
and N II 3995 A lines. Moreover, two O II lines, namely 
4093 A and 4705 A, were excluded because they show also 
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Table 2. List of the used C II, N II and O II lines 



CII 


El, 


logs/ 


Nil 


E u 


log gf 


Oil 


Ei, 


logs/ 


A, A 


eV 




A, A 


eV 




A, A 


eV 




3918.969 


16.332 


-0.533 


4035.080 


23.125 


0.623 


4072.150 


25.650 


0.552 


3920.681 


16.333 


-0.232 


4041.310 


23.142 


0.853 


4075.859 


25.665 


0.692 


5137.253 


20.702 


-0.940 


4043.533 


23.132 


0.743 


4085.114 


25.650 


-0.188 


5139.175 


20.704 


-0.740 


4199.977 


23.246 


0.030 


4132.800 


25.832 


-0.066 


5143.495 


20.704 


-0.212 


4227.738 


21.600 


-0.068 


4156.533 


25.849 


-0.696 


5145.165 


20.710 


0.189 


4447.028 


20.409 


0.285 


4185.451 


28.358 


0.604 


5151.083 


20.710 


-0.179 


4601.478 


18.466 


-0.428 


4303.821 


28.822 


0.640 


5648.070 


20.704 


-0.424 


4607.154 


18.462 


-0.507 


4317.134 


22.966 


-0.385 


5662.456 


20.710 


-0.249 


4613.870 


18.466 


-0.665 


4395.939 


26.249 


-0.169 


6578.048 


14.449 


-0.026 


4621.394 


18.466 


-0.514 


4414.889 


23.442 


0.171 


6582.873 


14.449 


-0.328 


4630.535 


18.483 


0.094 


4452.342 


23.442 


-0.789 


6779.940 


20.704 


0.024 


4643.083 


18.483 


-0.359 


4590.978 


25.661 


0.350 


6783.904 


20.710 


0.304 


4678.136 


23.572 


0.434 


4610.197 


29.063 


-0.170 








4987.377 


20.940 


-0.555 


4638.854 


22.966 


-0.332 








5002.700 


18.462 


-1.021 


4641.813 


22.980 


0.054 








5005.147 


20.666 


0.592 


4649.139 


22.999 


0.308 








5007.330 


20.940 


0.171 


4650.839 


22.966 


-0.361 








5010.619 


18.466 


-0.606 


4661.632 


22.980 


-0.277 








5025.654 


20.666 


-0.546 


4673.735 


22.980 


-1.088 








5045.092 


18.483 


-0.407 


4676.237 


22.999 


-0.395 








5666.628 


18.466 


-0.045 


4701.180 


28.830 


0.088 








5679.554 


18.483 


0.250 


4703.163 


28.513 


0.262 








5686.213 


18.466 


-0.549 


4890.865 


26.305 


-0.436 








6482.043 


18.497 


-0.162 


4906.841 
4941.064 
4943.000 


26.305 
26.554 
26.561 


-0.160 
-0.054 
0.239 




14000 16000 18000 20000 22000 24000 26000 28000 30000 32000 



Figure 4. Difference between the N abundance derived from the 
N II 3995 A line and the mean N abundance obtained from all 
other N II lines. 



a systematic trend with T e ff . It is interesting that in the re- 
cent work of Nieva & Przybilla (2012), where updated model 
atoms for carbon and nitrogen were applied, the lines C II 
4267 A and N II 3995 A did not show a marked difference 
from other lines. Therefore, the above-mentioned discrepan- 
cies for the lines C II 4267 A and N II 3995 A and two O II 
lines are likely to be due to some incompleteness of the C, 
N and O model atoms used in our analysis (Korotin et al. 
1998, 1999 a,b). 



5 EFFECTIVE TEMPERATURES T off AND 
SURFACE GRAVITIES log g 

Effective temperatures T c h and surface gravities logg of 
more than 100 B stars have been determined by us in Paper 
II. Now, we redetermine these parameters for the selected 
22 stars using the new stellar parallaxes of van Leeuwen 
(2007) as principal indicators of logg. Other indicators of 
log<? used in Paper II, namely the /3-index and equivalent 
widths of the Balmer lines H,3 and H 7 , are considered now 
only for comparison. As far as the effective temperature T e g 
is concerned, we use the same two indicators as in Paper II, 
i.e., the colour indices Q and [ci] in the photometric sys- 
tems UBV and uvby, respectively (as known, both indices 
are insensitive to interstellar extinction). 

Observed Q values were determined with photometric 
data from Mermilliod & Mermilliod's (1994) and Mermil- 
liod's (1994) catalogues. Observed [ci] and /? indices are 
found from Hauck & Mermilliod's (1998) catalogue. Ob- 
served equivalent widths W^H^) and W (H 7 ) are provided in 
Paper I. These observed parameters are compared with the 
computed values. Computations of W(H^) and W(H 7 ) for a 
large number of model atmospheres have been presented by 
Kurucz (1993). Colour indices for calculations of Q and [ci] 
are taken by us from Castelli & Kurucz (2003). Computed 
/3 values are published by Castelli & Kurucz (2006). 

As a typical example of the combined T e ff and loggr 
determination, we show in Fig. 5 the diagram for the star 
HR 1810. This is the standard T e ff - log g diagram, where 
various curves are the loci from the theoretical models which 
predict the observed value for each of the considered pa- 
rameters (Q, [ci], etc.). The locus that is set by parallax n 
is virtually a straight line; it depends on T e s very slightly. 
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Figure 5. The T eB - log g diagram for the star HR 1810. The filled 
circle corresponds to the adopted parameters T c ff = 21040 K and 
logg = 4.11. 
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Note that application of ir was described in detail by us ear- 
lier (Lyubimkov et al. 2009, 2010). In the case of HR 1810 
the parallax is rather large and known with a good accuracy, 
namely it — 5.20±0.21 mas. When the effective temperature 
T e ff is fixed, the surface gravity of the star is derived reliably 
from 7r; in fact, the uncertainty in logg is only ±0.03. 

Two indicators of T e g, namely indices Q and [ci], show 
in Fig. 5 somewhat different temperatures T c g; the difference 
is about 1000 K. We adopted the mean T c g value, so the final 
parameters of the star HR 1810 are T cff = 21040±520 K 
and logg = 4.11±0.07 (filled circle in Fig. 5). Fig. 5 shows 
that there is good agreement between the value logg = 4.11 
derived from the parallax and the logg estimates from the 
Balmer lines and /3-index. Using this technique, we obtained 
the parameters T e ff and log g for all programme stars; these 
values and their uncertainties are presented in Table 1. 

Fig. 5 displays one feature in the T e g determination 
for early B stars, which has been already noted by us in Pa- 
per II, namely: there is a systematic discrepancy between in- 
dices Q and [ci] as indicators of the T" e ff ■ An especially hotter 
Teff from Q as compared with [ci] is found for B stars with 
T ff > 24000 K, i.e., for spectral subtypes B0-B1.5. Our con- 
sideration of this discrepancy for B0-B1.5 stars in Paper II 
led to the conclusion that the T e ff([ci]) values are more cor- 
rect than T ff(Q). Theoretical Q values based on Kurucz's 
(1993) computed indices (U — B) and (B-V) have been cor- 
rected in order to eliminate disagreement between T c a(Q) 
and T e ff([ci]). Therefore, the index [ci] has been adopted as 
a principal indicator of T e ff. Now that we use the computed 
(U-B) and (B-V) indices from Castelli & Kurucz (2003), 
we find that the discrepancy between T c a(Q) and T e ft ([ci]) 
decreases but does not disappear. 

One may see from Table 1 that our programme B stars 
have effective temperatures from 15360 to 24090 K. For these 
relatively cool B stars the difference between T c g(Q) and 
T e ff ([ci]) leads to an uncertainty in T e ff of between ±20 and 
±870 K with a mean value of ±500 K. When comparing 
such accuracy with data on T c ff from other authors, one 
may conclude that the accuracy is plenty good enough. 

It is interesting to compare the T e g and logg values 




3.4 3.5 3.6 3.7 3.8 3.9 4.0 4.1 4.2 4.3 

logg (new values) 

Figure 6. Comparison of the new Teff and logg values with old 
parameters from Paper II. 

from Table 1 with the old T e ff and log g values for the same 
22 stars from Paper II. This comparison is shown in Fig. 6. 
One sees that new effective temperatures T e ff are slightly 
hotter than the old ones (upper panel). The mean difference 
is about 3 per cent, but for three hottest programme stars 
with T e ff close to 24000 K it reaches 4-5 per cent. This T e ff 
increment is mainly explained by the fact that now we added 
the second independent T e ff indicator, namely the index Q, 
which tends to increase T e ff in comparison with [ci]. 

Lower panel of Fig. 6 presents the comparison of new 
and old logg values. It is important to emphasize that 
our new surface gravities are found from stellar parallaxes 
and, therefore, are independent of model atmospheres of the 
stars. On the contrary, when the old logg values in Paper II 
were determined, the /3-index and equivalent widths of the 
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Figure 7. Comparison of our T e g values with ones derived with 
Napiwotzki's (2004) code. The straight line corresponds to equal 
effective temperatures. 



Balmer lines and H 7 were used; their observed values 
were compared with results of model computations. In other 
words, the new and old \ogg values are derived by different 
methods. One may see that there is a marked difference 
between the new and old log g values for some stars. For in- 
stance, the stars HR 1640 and HR 8385 (marked in Fig. 6) 
show a difference of +0.30 and -0.26 dex, respectively. Errors 
in their new log(? values are ±0.15 and ±0.10 dex, respec- 
tively (Table 1), so they can explain the discrepancy only 
partially. We believe that the new surface gravities log g de- 
rived from parallaxes are more reliable. It should be noted 
that the mean error of the logg values in Table 1 is ±0.11 
dex, whereas an accuracy of ±0.2 or ±0.3 dex was consid- 
ered typical for B stars a few years ago. Thus, our surface 
gravities logg derived from stellar parallaxes are rather ac- 
curate. 

What may we say about the accuracy of the derived ef- 
fective temperatures T e ff? To answer this question, we com- 
pare our T e g values with evaluations based on two other 
accurate T e s scales. Napiwotzki et al. (1993) provided the 
T c ff calibration for [u — b], (b — y)o and (B — V)o using a 
sample of stars with well known temperatures. Later Napi- 
wotzki (2004) presented a new calibration, which is based on 
colors derived from Kurucz's (1993) model atmospheres. A 
scale provided by Fitzpatrick & Massa (2005), which gives 
T c ff values in very good agreement with those derived with 
Napiwotzki's (2004) code, is based on a new calibration of 
synthetic photometry in various photometric systems for 45 
normal nearby B and early A stars on the basis of Kurucz's 
(1993) model atmospheres. Fitzpatrick & Massa emphasize 
that their and Napiwotzki's calibrations 'have no input data 
in common'. In fact, these calibrations use different observ- 
ables: Napiwotzki (2004) applies uvby/3 photometry, whereas 
Fitzpatrick & Massa (2005) employ IUE spectrophotometry, 
V-band fluxes and Hipparcos parallaxes (old catalog, ESA 
1997). 



Table 3. Comparison of our parameters with those of Nieva 
& Simon-Diaz (2011) = NST1 for two common stars, namely 
HR 1781 and HR 1820 (the NS'll loge(O) values are taken from 
Simon-Di'az, 2010) 



Parameter HR 1781 = HD 35299 
Present work NS'll 



HR 1820 = HD 35912 
Present work NS'll 



T eff , K 


23700±410 


24000±200 


19370±380 


19000±300 


logs 


4.27±0.07 


4.20±0.08 


4.05±0.19 


4.00±0.10 


14, km s~ x 


1.2±1.0 


Oil 


0.8± 1.0 


2±1 


loge(C) 


8.25±0.10 


8.37±0.07 


8.28±0.10 


8.33±0.09 


loge(N) 


7.76±0.08 


7.81±0.07 


7.76±0.09 


7.76±0.07 


loge(O) 


8.68±0.07 


8.71±0.08 


8.71±0.08 


8.77±0.12 



Table 4. Comparison of our parameters and the C, N and O 
abundances with those of Nieva & Przybilla (2012) = NPT2 for 
four common stars 



Parameter 


HR 1781 = 


HD 35299 


HR 1810 = 


HD 35708 




Present work 


NPT2 


Present work 


NP'12 


Teff, K 


23700±410 


23500±300 


21040±520 


20700±200 


log g 


4.27±0.07 


4.20±0.05 


4.11±0.07 


4.15±0.07 


Vt, km s" 1 


1.2±1.0 


Oil 


0.0±1.0 


2±1 


Ioge(C) 


8.25±0.10 


8.35±0.09 


8.25±0.12 


8.30±0.09 


loge(N) 


7.76±0.08 


7.82±0.08 


8.16±0.07 


8.22±0.07 


Ioge(O) 


8.68±0.07 


8.84±0.09 


8.65±0.06 


8.82±0.11 




HR 2928 = 


HD 61068 


HR 8385 = 


HD 209008 




Present work 


NPT2 


Present work 


NP'12 


Teff, K 


23650±850 


26300±300 


15360±240 


15800±200 


logg 


3.83±0.11 


4.15±0.05 


3.76±0.10 


3.75±0.05 


Vt, km s- 1 


3.7±1.0 


3±1 


1.2±1.0 


4±1 


loge(C) 


8.01±0.10 


8.27±0.07 


8.57±0.19 


8.33±0.09 


loge(N) 


7.91±0.10 


8.00±0.12 


7.97±0.12 


7.80±0.11 


Ioge(O) 


8.69±0.11 


8.76±0.09 


8.92±0.12 


8.80±0.11 



We show in Fig. 7 a comparison of our T c ff values from 
Table 1 with ones derived by us from wby/3-photometry 
(Hauck & Mermilliod 1998) with Napiwotzki's (2004) code. 
One may see that there is very good agreement with Napi- 
wotzki's T e ff scale and, therefore, with Fitzpatrick & Massa's 
scale. Note especially that no systematic discrepancy is seen. 
This agreement with these two independent accurate T e s 
scales confirms the reliability of our new T c ff values. 

Additional checks on the parameters are possible 
through comparisons for stars common to recent studies. 
For example, Nieva & Simon-Di'az (2011) studied 13 early 
B-type stars in the Ori OBI association including deriva- 
tions of the C, N and O abundances. Their sample includes 
two stars from Table 1, namely HR 1781 and HR 1820. We 
compare in Table 3 the parameters T ef j and logg of these 
stars from their and our work. One may see that there is very 
good agreement with the T e g and log g values; the differences 
are within errors of the T e g and logg determinations. It is 
important to note that Nieva & Simon-Di'az used quite dif- 
ferent methods to us for their T c s and logg determinations, 
i.e., their methods are based on the non-LTE analysis of 
lines of C, O, Si and Fe in two stages of ionization. So, the 
comparison with their data confirms the reliability of our 
new T e ff and log g values. Nieva & Simon-Di'az adopt as we 
do the Kurucz ATLAS model atmospheres. 

Nieva & Przybilla's (2012) study of 29 early B stars 
contains four objects in common with our list. (Note that 
this 2012 paper determines T c ff and logg by the methods 
very similar to those used by Nieva & Simon-Di'az, 2011). 




Figure 8. Claret's (2004) evolutionary tracks and positions of 
programme stars on the log T e g - log g plane. Error bars in log T c g 
and logg are shown for each star. 



According to Table 4, there is good agreement between the 
parameters Teff and logg of these stars and our values, ex- 
cept for HR 2928; for this star Nieva & Przybilla (2012) 
obtained a markedly hotter effective temperature T e e and, 
as a result, the greater \ogg value. It should be noted that 
Napiwotzki's (2004) method gives for HR 2928 the parame- 
ters T ft = 23980 K and logg = 3.83 that are very close to 
ours. 

When basic parameters T e ff and logg are known, one 
may determine the mass M of a star from evolutionary 
tracks. The derived masses M of the programme stars (in the 
solar masses M©) are presented in Table 1; they are found 
from Claret's (2004) evolutionary tracks. One sees that the 
masses range from 5.4 to 11.2 Mq. The typical error in the 
M determination due to uncertainties in T e g and log g is ±8 
per cent. 

A question may arise about the reliability of stellar 
masses M ev derived from evolutionary tracks. It is known 
that the most reliable masses M or t are obtained for compo- 
nents of binaries from an analysis of their orbital elements. 
Lyubimkov (1996) compared the M ev and M or b values for 
a wide mass range; one may see from this comparison that 
there is very good agreement between M ev and M or t, for 
masses between 5 and 15 Mq (see Fig. 9 there). So, we be- 
lieve that our M values in the 5-11 Mq range inferred from 
Claret's (2004) evolutionary tracks are rather accurate. 

In Fig. 8 we show Claret's (2004) tracks in the T e g-\ogg 
plane, as well as locations of all 22 stars with corresponding 
error bars for T e a and logg; note that the labels for the 
tracks are log M /Mq . One sees that all the stars are inside 
the MS band. Two stars, namely HR 7426 and HR 7996 with 
the lowest surface gravity logg = 3.60, are close to the MS 
phase termination. Nonetheless, as mentioned in Section 2, 
the observed C, N and O abundances for such stars (they 
have M ~ 6-7 Mq and usini ~ 30 km s _1 ) are expected 
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Table 5. Carbon, nitrogen and oxygen abundances in atmo- 
spheres of 22 B-type MS stars. The number of lines used is shown 
in brackets. 



HR 


V t 

km s^ 1 


loge(C) 




loge(N) 




loge(O) 




38 


2.7 


8.29±0.09 


(9) 


7.70±0.06 


(11) 


8.56±0.12 


(8) 


1617 


0.5* 


8.16±0.06 


(6) 


7.63±0.11 


(11) 


8.63±0.15 


(11) 


1640 


0.0 


8.36±0.13 


(4) 


7.78±0.13 


(11) 


8.83±0.11 


(9) 


1781 


1.2 


8.25±0.10 


(11) 


7.76±0.08 


(23) 


8.68±0.07 


(16) 


1810 


0.0 


8.25±0.12 


(9) 


8.16±0.07 


(13) 


8.65±0.06 


(8) 


1820 


0.8 


8.28±0.10 


(6) 


7.76±0.09 


(19) 


8.71±0.08 


(12) 


1923 


1.1 


8.25±0.06 


(9) 


7.74±0.08 


(20) 


8.62±0.06 


(12) 


1933 


1.8 


8.34±0.14 


(4) 


7.81±0.12 


(12) 


8.58±0.10 


(6) 


2058 


0.8 


8.25±0.09 


(9) 


7.72±0.10 


(22) 


8.68±0.09 


(16) 


2205 


1.9 


8.33±0.09 


(8) 


7.76±0.07 


(20) 


8.66±0.07 


(17) 


2344 


0.0* 


8.21±0.02 


(3) 


7.88±0.20 


(5) 


8.73±0.08 


(5) 


2756 


0.0* 


8.46±0.19 


(5) 


7.82±0.03 


(1) 


8.91±0.07 


(4) 


2824 


1.0 


8.23±0.10 


(8) 


7.66±0.09 


(22) 


8.63±0.12 


(20) 


2928 


3.7 


8.01±0.10 


(7) 


7.91±0.10 


(19) 


8.69±0.11 


(16) 


3023 


1.2 


8.30±0.12 


(9) 


7.79±0.13 


(16) 


8.68±0.11 


(12) 


7426 


0.0* 


8.30±0.16 


(9) 


7.65±0.13 


(13) 


8.59±0.14 


(9) 


7862 


2.0* 


8.48±0.12 


(3) 


7.86±0.08 


(6) 


8.95±0.10 


(2) 


7996 


0.1* 


8.53±0.14 


(3) 


7.87±0.09 


(7) 


8.84±0.09 


(4) 


8385 


1.2* 


8.57±0.19 


(9) 


7.97±0.12 


(11) 


8.92±0.12 


(8) 


8549 


1.6 


8.26±0.07 


(8) 


7.74±0.07 


(20) 


8.70±0.07 


(23) 


8768 


1.4 


8.42±0.10 


(12) 


7.86±0.07 


(20) 


8.88±0.16 


(18) 


9005 


2.1 


8.29±0.09 


(10) 


7.84±0.08 


(22) 


8.59±0.07 


(21) 



*) These Vt values are determined in Paper III from He I lines 



to be close to the initial ones; rotationally-induced mixing 
is not predicted to affect surface abundances unless a slow 
rotator is a rapid rotator seen at a high angle of inclination. 



6 CARBON, NITROGEN AND OXYGEN 
ABUNDANCES 

As mentioned above, model atmospheres for the stars were 
computed using Kurucz's (1993) code ATLAS 9 for the de- 
rived parameters T e ff and logg. When determining element 
abundances, apart from the star's model atmosphere the 
microturbulent parameter Vt should be determined. We in- 
ferred the Vt values from N II and O II lines. A standard 
procedure was used, namely: we searched for the Vt that 
yields a zero slope for the abundance vs. W relation; in other 
words, no trend of the N and O abundances with equivalent 
widths W must exist. The derived parameters Vt (N II) and 
Vt (O II) were applied to the N and O abundance analy- 
ses, respectively. Since the few available C II lines are weak, 
the C abundance is only very slightly dependent on Vt; we 
adopted the mean of the Vt (N II) and Vt (O II) - see Ta- 
ble 5. For the coolest programme stars, a reliable derivation 
of Vt from N II and O II lines is impossible, because the few 
available lines are weak. Therefore, for these stars (seven in 
all) we have taken Vt from He I lines (see Paper III); these Vt 
values are marked in Table 5 by asterisks. Table 5 shows that 
the Vt values for the programme stars are rather small: they 
vary from 0.0 to 3.7 km s _1 with the average of 1.1 km s _1 . 
Of course, when all lines are weak, the derived abundance 
is independent of the adopted value of Vt. 

The derived C, N and O abundances with the standard 
deviations from the line-to-line scatter for 22 programme 
stars are presented in Table 5. The number of C II, N II and 
O II lines used for each star is given there (in brackets). This 
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number depends on the effective temperature T c g, and rota- 
tional velocity vsini, as well as the quality of the observed 
spectrum of a star. The C, N and O abundances are given in 
Table 5 as the mean values log t with the standard deriva- 
tion a. On average a = ±0.11 dex for C and ±0.10 dex for 
N and O. Averaging data from Table 5, we obtain the fol- 
lowing mean C, N and O abundances: loge(C) = 8.31±0.13, 
loge(N) = 7.80±0.12 and loge(O) = 8.73±0.13. 

As mentioned above, the accuracy of the C, N and 
O abundances depends on the reliability of the basic pa- 
rameters T c ff and logg. In order to estimate the accuracy, 
one may compare the derived mean loge(C), loge(N) and 
log e(O) values with results obtained only for the stars with 
the most reliable parameters T e g and log jr. First, there 
are 11 stars in Table 1 with the relatively small tem- 
perature errors AT c ff < 500 K; we found for them on 
average loge(C) = 8.34±0.13, loge(N) = 7.78±0.09 and 
loge(O) = 8.74±0.11. Second, the most accurate loggr val- 
ues are obtained for the stars with small errors in the par- 
allaxes 7r. There are 7 programme stars with the mini- 
mum errors in tt, i.e., errors of 10 per cent. The mean 
abundances for these 7 stars are loge(C) = 8.28±0.10, 
loge(N) = 7.81±0.18 and loge(O) = 8.69±0.10. Third, there 
are 4 common stars from both previous groups, which have 
simultaneously the most accurate T e g and log 3 values; we 
found for them loge(C) = 8.30±0.13, loge(N) = 7.74±0.08 
and loge(O) = 8.72±0.13. One sees that in all three cases 
the mean C, N and O abundances virtually coincide with the 
above-mentioned averages for a total sample of 22 stars. It 
shoud be noted that van Leeuwen (2007) did not recommend 
use of parallaxes n with errors > 10 per cent. However, we 
found no marked difference in the mean C, N and O abun- 
dances between objects with errors in n > 10 per cent (15 
stars) and ^ 10 per cent (7 stars): the difference is 0.04, 0.00 
and 0.05 dex for C, N and O respectively, and these values 
are less than uncertainties in the mean abundances. 

We estimated typical errors in the derived C, N and O 
abundances, which originate from four contributors, namely 
uncertainties in the parameters T e ff, log<? and Vt and the 
scatter in the log e estimates from various lines. Since equiv- 
alent widths of C II, N II and O II lines depend nonlinearly 
on Teff (see, e.g., Figs 2 and 3), the contribution of the un- 
certainty in T e ff to a total error can be dependent on T e g. 
We selected a number of the stars from Table 1, which rep- 
resent the whole T c g range, and estimated for them the total 
errors in loge(C), loge(N) and loge(O) adopting the follow- 
ing mean uncertainties: AT e ff = ±500 K, A log g = ±0.11 
and AVt = ±1.0 km s _1 ; the scatter in loge between lines 
is characterized by the mean dispersion a = ±0.10 dex (see 
above). We found that the total error in log e(C) varies from 
0.17 to 0.11 dex when T eS varies from 15000 to 24000 K; 
the total error in loge(N) and loge(O) varies from 0.18 to 
0.11 dex and from 0.21 to 0.13 dex, respectively. The con- 
tribution of the uncertainty AVt = ±1.0 km s _1 seems to 
be imperceptible for all temperatures T e g as the lines used 
are rather weak. The contribution of the T e g uncertainty 
increases with decreasing T e g. For the hottest stars of our 
sample (T c g ~ 22000-24000 K) the errors in log e originate 
mostly from the scatter in loge between lines, whereas for 
the cooler stars (T cff ~ 15000-19000 K) the contribution of 
the T e ff uncertainty is significant. 

Independent but limited confirmation of the accuracy 



Table 6. Comparison of our mean C, N and O abundances with 
other recent results 



Source 


Number 
of stars 


loge(C) 


loge(N) 


loge(O) 


Present work 


22 


8.31±0.13 


7.80±0.12 


3.73±0.13 


Takeda ct al. (2010) 


64 






3.71±0.06 


Nieva & 


13 a > 


8.35±0.03 


7.82±0.07 


3.77±0.03 


Simon-Diaz (2011) 










Cunha, Hubeny & 


10") 






3.78±0.05 


Lanz (2012) 










Nieva & Przybilla (2012) 


29 


8.33±0.04 


7.79±0.04 


3.76±0.05 



a ) stars in the Ori OBI association 

of our C, N and O abundances is presented in Table 3, 
where our and Nieva & Simon-Diaz's (2011, hereinafter 
NST1) data for two common stars from the Ori OBI as- 
sociation are compared. The derived C, N and O abun- 
dances show very good agreement: differences in the log e(C), 
log e(N) and log e(O) values vary from 0.00 to 0.06 dex, i.e. 
they are markedly less than errors of the loge(C), loge(N) 
and loge(O) determination with only the C abundance for 
HR 1781 showing a greater difference, namely 0.12 dex. 

In the same vein, we compare in Table 4 our C, N and 
O abundances with data of Nieva & Przybilla (2012, here- 
inafter NPT2) for four common stars. In particular, one sees 
from Table 4 that the above-mentioned enhanced N abun- 
dance for HR 1810 is confirmed by NPT2. In general, an 
agreement between our and NPT2 results is worse than for 
two stars from NST1 (Table 3). Especially large discrepan- 
cies with NPT2 are found for the C abundances in HR 2928 
and HR 8385 (±0.26 and -0.24 dex, respectively). Note that 
in the case of HR 2928 the discrepancy is partially explained 
by the substantial difference between the adopted T e ff val- 
ues. 

It should be noted that the C, N and O abundances 
derived in NST1 and NPT2 show no correlation with the 
effective temperature T e ft (see Fig. 7 in NPT2). Our N and 
O abundances show no dependence on T e g, too. However, 
unlike N and O, the C abundances for cooler stars of 
our sample with T e ft ^ 18000 K seem to be somewhat 
higher than for hotter stars. In fact, the mean value is 
loge(C) = 8.46±0.09 for the stars with T cB < 18100 K 
(6 objects) and loge(C) = 8.25±0.08 for the stars with 
T c ff > 18500 K (16 objects). Note that there are only two rel- 
atively cool B stars in NS'll and NPT2 with T cff < 18000 K; 
their C abundances derived from C II lines are in good agree- 
ment with the loge(C) values for other stars. On the one 
hand, marked difference in the C abundances between the 
stars of our sample with T cff < 18100 K and T cff > 18500 K 
can result from some incompleteness of the used C model 
atom (in NS'll and NPT2 the updated model is applied). 
On the other hand, it is interesting that the mean value 
log e(C) = 8.46 for the cooler stars agrees very well with the 
solar C abundance, whereas the mean value log e(C) = 8.25 
for the hotter stars (as well as the mean C abundances in 
NS'll and NPT2) show a marked underabundance. 



7 DISCUSSION 

When our work was in progress, a several interesting publi- 
cations appeared with C, N and O abundances in B stars. In 



Carbon, Nitrogen and Oxygen Abundances 11 



Table 7. The mean C, N and O abundances for the B-type MS 
stars in comparison with the solar abundances 



Object loge(C) loge(N) loge(O) Source 

B-typc stars 8.31±0.13 7.80±0.12 8.73±0.13 Present work 

Solar 8.50±0.06 7.86±0.12 8.76±0.07 Caffau et al. 

photosphere (2008, 2009, 2010) 

Solar 8.43±0.05 7.83±0.05 8.69±0.05 Asplund ct al. 

photosphere (2009) 

Protosolar 8.47±0.05 7.87±0.05 8.73±0.05 Asplund ct al. 

abundances (2009) 



Table 6, we compare our mean C, N and O abundances with 
these recent results. Some of the 64 B stars from Takeda et 
al.'s (2010) sample are close to or beyond the termination 
of the MS (see their Fig.l). However, due to the relatively 
low masses (M = 3-9 Mq) and slow rotational velocities 
(vsini ^ 30 km s -1 ), these stars are expected to have re- 
tained their initial surface C, N and O abundances. Takeda 
et al. determined the O abundance from the O I 6156 A 
multiplet, and the Ne abundance from Ne I lines but did 
not provide either C or N abundances. Of the 29 early B 
stars analysed by NP'12, most are caught well before the 
end of the MS (see Fig. 5 there). An obvious exception is 
the star a Pyx (HR 3468, HD 74575) with the lowest surface 
gravity log(? = 3.60, which is close to the MS's termination. 
Its mass is M ~ 12 Mq and the projected rotational veloc- 
ity is vsini = 11 km s _1 ; according to (Frischknecht et al. 
2010) and assuming the star is not a rapid rotator observed 
at a high angle of inclination, such a star cannot change 
markedly its surface C, N and O abundances during the MS 
phase. On the whole, one may suppose that the vast ma- 
jority of the C, N and O abundances presented in Table 6 
correspond to the initial abundances in these young stars. 

Table 6 shows that our mean C, N and O abundances 
are in very good agreement with data from recent papers. 
In particular, there is excellent agreement with the mean 
C, N and O abundances found in NS'll and NP'12. This is 
very important because both in the T e g and logg determi- 
nation and in the CNO abundance analysis quite different 
techniques were used by NS'll (and NPT2) and us. 

It should be noted that NS'll, as well as Cunha, 
Hubeny & Lanz (2012), studied the unevolved early B stars 
in the Ori OBI association. In both studies, the T e g range 
of the stars is displaced to hotter temperatures relative to 
our sample: 19000-33400 K for NS'll against 15360-24090 K 
for our sample. The agreement over abundances extends to 
Mg: 7.57±0.06 from NS'll and 7.59±0.15 from Paper IV. 

Since the Ori OB 1 stars have been formed from in- 
terstellar material of the (presumed) same metallicity, it is 
quite natural that they show very homogeneous initial C, 
N and O abundances. In contrast, the scatter in the indi- 
vidual loge(C), loge(N) and loge(O) values for the stars of 
our sample is greater than for the above samples from Ori 
OB 1. However, the stars of our sample do not belong to a 
single cluster or association; on the contrary, they are dis- 
tributed around the solar neighbourhood up to d — 600 pc. 
Therefore, we speculate that the spread among the individ- 
ual loge(C), loge(N) and loge(O) values for our stars may 
be augmented by real star-to-star variation in the initial C, 
N and O abundances; the Galactic abundance gradient may 
be one source of the variation. 




0.0 r 1 

0.0 0.1 0.2 0.3 0.4 

N/O 

Figure 9. Relation between the mass ratios N/C and N/O for 22 
programme stars. Two stars with the high N/C values, namely 
HR 1810 and HR 2938, are marked. Theoretical path of the N/C 
vs. N/O relation is shown by the dashed line drawn through the 
initial values N/O = 0.103 and N/C = 0.342 (open star in the 
figure) with a slope d(N/0)/d(N/C) = 4.6 from NP'12. 

In Table 7 our mean log e(C), log e(N) and log e(O) val- 
ues for 22 B stars are compared with the recent solar val- 
ues (Caffau et al. 2008, 2009, 2010; Asplund et al. 2009). 
Protosolar abundances are 0.04 dex larger (Asplund et al. 
2009). When comparing the results in Table 7, we conclude 
that our mean log e(N) and log e(O) values for the B stars 
coincide within errors of their determination with the pho- 
tospheric and protosolar abundances. However, our mean 
loge(C) value, also those in NS'll and NP'12 (Table 6), 
shows a marked underabundance in comparison with the 
Sun (the deficiency is especially noticeable when compared 
with small errors in the mean C abundances in NS'll and 
NP'12). 

We also determined the ratios N/O and N/C for 22 
programme stars; we found the mean values log (N/O) = - 
0.91±0.13 and log (N/C) = -0.50±0.15. The first value is 
very close to the solar ratio log (N/O) = -0.90±0.14 (Caf- 
fau et al.) or -0.86±0.07 (Asplund et al.). The second one 
is somewhat greater than the solar ratio log (N/C) = - 
0.64±0.13 (Caffau et al.) or -0.60±0.07 (Asplund et al.); 
a main cause of this difference is a slight C deficiency in B 
stars. 

Following to NP'12 (see their Fig. 14), we show in Fig. 9 
the N/C vs. N/O diagram, where the derived CNO abun- 
dances, like NP'12, have been converted into the mass frac- 
tion scale. One may see that there is no correlation between 
N/O and N/C for the most of the stars; this is explain- 
able because their observed CNO abundances are expected 
to correspond to the unchanged initial values. However, two 
of the stars, namely HR 1810 and HR 2938, show the rel- 
atively high N/C values. One may suppose that these two 
stars could change their atmospheric abundances of light ele- 
ments during the MS phase evolution due to the rotationally 
induced mixing. However, their positions in the T e a - log<? 
diagram (Fig. 8) are not exclusive. Moreover, the presented 
errors in the N/O and N/C values for these stars are large 
(note that the errors in Fig.9 are proportional to N/O and 
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Table 8. Comparison of our mean C, N and O abundances for 20 
B stars (the 'mixed' stars HR 1810 and HR 2928 are omitted) with 
a present-day Cosmic Abundance Standard (CAS) from NP'12 



Element 


Our abundances 


CAS 






(NP'12) 


C 


8.33±0.11 


8.33±0.04 


N 


7.78±0.09 


7.79±0.04 


O 


8.72±0.12 


8.76±0.05 



N/C themselves), so their difference from other stars can- 
not recognize as a reliable one. Nevertheless, it should be 
noted that the relatively high N/C ratios for HR 1810 and 
HR 2938 were confirmed by NP'12. 

We also display in Fig. 9 a possible theoretical trend 
(dashed line) which corresponds to the NP'12 slope 
d(N/0)/d(N/C) = 4.6 and our initial values N/O = 0.103 
and N/C = 0.342 found for 20 stars without HR 1810 and 
HR 2938 (open star in Fig. 9); note that these initial val- 
ues are very close to the NP'12 ones. When excluding for 
the assumptive mixed stars HR 1810 and HR 2938, we 
obtain the following mean abundances for the remaining 
20 stars: loge(C) = 8.33±0.11, loge(N) = 7.78±0.09 and 
loge(O) = 8.72±0.12. We compare in Table 8 these values 
with a present-day Cosmic Abundance Standard (CAS) of 
NP'12; one sees that there is an excellent agreement. 

Thus, our findings, as well as up-to-date data of other 
authors show that the observed (unevolved) N and O abun- 
dances in B stars are very close to the solar ones. On the 
other hand, carbon shows a small stable deficiency. In fact, 
using the old carbon model atom of Sigut (1996), we found 
the mean abundance loge(C) = 8.31±0.13. Basing on the 
updated carbon model atom, NS'll and NP'12 obtained 
the very close values 8.35±0.03 and 8.33±0.04, respectively 
(note that errors are very small there). All these values are 
markedly lower than the current log e(C) evaluations for the 
Sun (Table 7) . It is important to note that many other chem- 
ical elements show abundances that are in an excellent agree- 
ment with the solar ones (see Introduction). Why does only 
carbon show such a unique difference from the Sun's abun- 
dance? 

It is possible that there are some specific details in non- 
LTE computations of C II lines that were ignored up to now. 
In particular, the C II ionization energy is E[ on — 24.383 eV 
that corresponds to the wavelength A = 508.5 A; therefore, 
the C II photoionization is controlled by the EUV radia- 
tion. For example, such a radiation was observed for the star 
e CMa (B2 II); it was found that the observed EUV flux at 
A < 504 A (He I continuum) is significantly greater than 
the computed one (Cassinelli et al. 1995). Later Gregorio 
at al. (2002), using both the ATLAS9 model atmospheres 
of Kurucz (1993) and the spherical model atmospheres of 
Aufdenberg et al. (1999), showed that a more satisfactory 
agreement with the observed EUV flux for the star e CMa 
can be achieved; however, the observed flux has a different 
shape than the predicted one. Moreover, the observed line 
intensities in this region seem to be much stronger than pre- 
dicted ones. Thus, the theory cannot describe adequately 
the EUV flux which is important for computations of the 
C II photoionization (as well as the N II and O II pho- 
toionization, too). Moreover, we cannot exclude that some 



further update of the carbon model atom is needed to solve 
the problem. 

Another alternative was proposed recently; it was sup- 
posed that the Sun can migrate during its life (i.e., for 
the past 4.5 ■ 10 9 yr) from inner parts of the Galactic disk 
where it has born, so its observed chemical composition may 
differ from the composition of young stars in its present 
neighbourhood (typical ages of the B stars in question are 
t ~ (10 — 100) ■ 10 6 yr). This hypothesis is discussed, e.g., in 
NP'12 (see also references there). Detailed analysis of the ra- 
dial migration of stars (including the Sun) within the Galaxy 
is very hard; nevertheless, the question remains: why is it 
that only the carbon abundance shows a marked difference 
from the Sun's abundances? In our opinion, the slight carbon 
deficiency in the unevolved B stars is a subject for further 
discussion. 



8 CONCLUDING REMARKS 

We selected from our original list of early and medium B- 
type stars (Paper III) the relatively cool stars with low ro- 
tational velocities (vsini < 70 km s _1 ). Their two basic 
parameters, namely the effective temperature T e g and the 
surface gravity \ogg were redetermined; application of new 
stellar parallaxes allowed increasing markedly the accuracy 
of the derived logg values. Our final list contains 22 stars 
with the following parameters: (i) effective temperatures T e a 
are less than 24100 K; (ii) surface gravities logp are mostly 
greater than 3.75; (iii) the projected rotational velocities 
are vs'mi ^ 66 km s _1 . These stars with moderate masses 
M = 5-11 Mq and slow rotation should keep their initial 
surface C, N and O abundances according to the theory of 
rotationally-induced mixing. 

We found for these stars the mean abundances 
loge(C) = 8.31±0.13, loge(N) = 7.80±0.12 and 
loge(O) = 8.73±0.13; we suppose that these values 
are the initial CNO abundances in B-type MS stars. The 
derived N and O abundances coincide with the solar values, 
whereas the C abundance is somewhat lower than the solar 
one. We showed that these C, N and O abundances are in 
good agreement with the recent values for B-type MS stars 
by other authors. In particular, the small C deficiency is 
confirmed by the recent studies by NS'll and NP'12. 

Two programme stars, namely HR 1810 and HR 2938, 
can be mixed during the MS phase. If these stars are omit- 
ted, we obtain the following mean abundances for the re- 
maining 20 stars: log e(C) = 8.33±0.11, log e(N) = 7.78±0.09 
and loge(O) = 8.72±0.12; these values are in excellent 
agreement with a present-day Cosmic Abundance Standard 
(CAS) of Nieva & Przybilla (2012). 

Thus, on the one hand, the initial N and O abundances 
in nearby B-type stars, as well as the abundances of some 
other metals, confirm that young stars in the solar neigh- 
bourhood and the Sun have the same metallicity. In the 
absence of a migration of the Sun from its birthplace, one 
may conclude that during the Sun's life (i.e., for the past 
4.5 ■ 10 9 yr) the metallicity of the solar neighbourhood has 
not markedly changed; so, an intensive enrichment of the 
solar neighbourhood by metals occurred before the Sun's 
birth. On the other hand, carbon is an evident 'outlier' in 
this scenario, because the observed C abundance in B stars 
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shows a small deficiency in comparison with the Sun. This 
carbon deficiency demands an alternative hypothesis either 
the Sun has migrated from inner parts of the Galactic disk 
where it has born or present non-LTE computations over- 
look one or more physical effects in the carbon ions. The C 
underabundance in B stars needs further study. 
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